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Abstract 

Continuum states of the Dirac equation are calculated numerically for the electrostatic field 
generated by the charge distribution of an atomic nucleus. The behavior of the wave functions of 
an incoming electron with a given asymptotic momentum in the nuclear region is discussed in detail 
and the results are compared to different approximations used in the data analysis for quasielastic 
electron scattering off medium and highly charged nuclei. It is found that most of the approximations 
provide an accurate description of the electron wave functions in the range of electron energies above 
100 MeV typically used in experiments for quasielastic electron scattering off nuclei only near the 
center of the nucleus. It is therefore necessary that the properties of exact wave functions are 
investigated in detail in order to obtain reliable results in the data analysis of quasielastic (e, e'p) 
knockout reactions or inclusive quasielastic (e, e') scattering. Detailed arguments are given that 
the effective momentum approximation with a fitted potential parameter is a viable method for 
a simplified treatment of Coulomb corrections for certain kinematical regions used in experiments. 
Numerical calculations performed within the framework of the single particle shell model for nucleons 
lead to the conclusion that our results are incompatible with calculations performed about a decade 
ago, where exact electron wave functions were used in order to calculate Coulomb corrections in 
distorted wave Born approximation. A discussion of the exact solutions of the Dirac equation for 
free electrons in a Coulomb field generated by a point-like charge and some details relevant for the 
numerical calculations are given in the appendix. 
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1 Introduction 

Quasielastic (e, e'p) knockout reactions provide a powerful possibility to obtain information on the 
electromagnetic properties of nucleons embedded in the nuclear medium, since the transparency of the 
nucleus with respect to the electromagnetic probe makes it possible to explore the entire nuclear volume. 
Inclusive (e, e') scattering, where only the scattered electron is observed, provides information on a 
number of interesting nuclear properties like, e.g., the nuclear Fermi momentum pQ, high-momentum 
components in nuclear wave functions 0, modifications of nucleon form factors in the nuclear medium 
[3], the scaling properties of the quasielastic response allow to study the reaction mechanism [I], 
and extrapolation of the quasielastic response to infinite nucleon number A = oo provides us with 
a very valuable observable of infinite nuclear matter There is now considerable theoretical and 
experimental interest in extracting longitudinal and transverse structure functions as a function of 
energy loss for fixed three momentum transfer for a range of nuclei 0. In August 2005, the Thomas 
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Jefferson Nationaf Accelerator Facility (TJNAF) Proposal E01-016, entitled 'Precision measurement 
of longitudinal and transverse response functions of quasi-elastic electron scattering in the momentum 
transfer range 0.55 GeV < \q\ < 1.0 GeV was approved such that the experiments will be performed 
in the near future at the TJNAF using 4 He, 12 C, 56 Fe and 208 Pb as target nuclei. 

The plane wave Born approximation is no longer adequate for the calculation of scattering cross 
sections in the strong and long-range electrostatic field of highly charged nuclei, and it has become 
clear in recent years that the correct treatment of the Coulomb distortion of the electron wave function 
due to the electrostatic field of the nucleus is unavoidable if one aims at a consistent interpretation of 
experimental data. E.g., it is still unclear whether the Coulomb sum rule is violated in nuclei [""". 

Distorted wave Born approximation (DWBA) calculations with exact Dirac wave functions have 
been performed by Kim et al. |B| in the Ohio group and Udias et al. (HI E] for quasielastic scattering 
off heavy nuclei. However, these calculations are cumbersome and difficult to control by people who 
were not directly involved in the development of the respective programs. Early DWBA calculations 
for 12 C and 40 Ca were presented in jllj 

Various approximate treatments have been proposed in the past for the treatment of Coulomb 
distortions "J"!"] El EEH El El El El El > and there is an extensive literature on the so-called eikonal 
approximation |20 [ I21 | 122*1 123 | "2" H I25[ l26j . At lowest order, an expansion of the electron wave function 
in aZ, where a is the fine-structure constant and Z the charge number of the nucleus, leads to the well 
known effective momentum approximation (EMA) [2"""|, which plays an important role in experimental 
data analysis and which will be explained below. 

The effect of the charged nucleus on the electron wave function is twofold: Firstly, the (initial and 
final state) electron momentum kij is enhanced in the vicinity of the nucleus due to the attractive 
electrostatic potential, i.e., the wave length of the electron is becoming shorter near the nucleus. 
Secondly, the attractive potential of the nucleus leads to a focusing of the electron wave function in 
the nuclear region. Solutions for the Dirac equation for the scattering of electrons in the nuclear 
field can be obtained from a partial wave expansion, where the radial Dirac equation must be solved 
numerically for each partial wave. To avoid such a computational effort, an approximate treatment is 
often adopted, based on a high-energy expansion in inverse powers of the electron energy |12| I13 | "15] . 
The resulting expression for the distorted electron wave function is then expanded in powers of aZ as 
(we use units with h = c = 1 throughout) 

A = e ±i5 ^je^[l + gU(a,b,k',r) + g® (a,b,k' ,r) + ...}u T , (1) 

where the sign ± refers to the two scattering solutions with outgoing or incoming spherical waves, 
respectively, whereas the corresponding indices i, f are neglected for the sake of notational simplicity, 
u T is the plane wave spinor for the electron with a given helicity r and k is the asymptotic electron 
momentum. For the case of a uniform spherical charge distribution of radius R, the values of k', a, b 
and 5 1/2 are given by jXJSJ 

k' = k + -^, 5 1/2 = aZ^--log2kR)+b, 
aZ 3aZ 

with k' parallel to k. These values enter the first-order term according to 

g 0) = ar 2 + iar 2 k'f± ib[(k' x f) 2 + 2ik'r - s(k' x r)], (3) 

where the spin operator s describes spin-dependent effects which are comparably small for higher 
electron energies. The meaning of the parameter a can be easily understood from a semiclassical 
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observation. For a highly relativistic electron with mass m falling along the z-axis (i.e. with zero 
impact parameter and > V > m) towards the nuclear center, the momentum inside the spherical 
charge distribution is given by 

iW = *-„(„, V(r^-^ + §(£ (4) 

where V(r) is the electrostatic potential inside the charged sphere. Modifying the electron plane wave 
phase e lkr = e lkz to e tk ' z i l+az2 ) ; as it is induced at lowest order by the second term in leads to the 
z— dependent electron momentum k(z) 

ld ik , z{1+az2) = y + = ~ k{z yk>z(i+az*) (5) 

i dz 

in agreement with eq. @. The parameter b describes mainly the deformation of the wave front and 
can also be derived from semiclassical observations. 

The standard method (in the case of light nuclei) to handle Coulomb distortions for elastic scattering 
in data analysis is the effective momentum approximation (EMA), which corresponds to the lowest 
order description of the Coulomb distortion in aZ. EMA accounts for the two effects of the Coulomb 
distortion mentioned above (momentum modification and focusing) in the following way. For a highly 
relativistic electron with zero impact parameter the so-called effective momenta k[ ^ of the electron are 
given by 

k'i = ki + Ak, kf = kf + Ak, kij = \h,f\, Kj = Wi f\i Ak = —Vo/c, (6) 

where Vo is the potential energy of the electron in the center of the nucleus in analogy with eq. |2J). 
E.g., for 208 Pb we have Vo ~ — 25MeV, not a negligible quantity when compared to energies of some 
hundreds of MeV typically used in electron scattering experiments. Cross sections are then calculated 
by using plane electron waves corresponding to the effective momenta instead of the asymptotic values 
in the matrix elements, and additionally one accounts for the focusing factors of the incoming and 
outgoing electron wave k[/ki and k'f/kf, which both enter quadratically into the cross sections. The 
main problem of the method is the fact that both the focusing and the electron momentum are not 
constant inside the nuclear volume. In the case of nucleon knockout reactions, most of the hit nucleons 
are located near the surface of the nucleus, where the classical momentum of the electrons is not given 
by the central value. 

A strategy to remedy this defect is to alter the definition of the effective momenta by not using 
the central potential value Vo, but a value V(ry) obtained from some fitting procedure (see also |28| 
and references therein). An even more ambitious strategy would be the introduction of two effective 
momenta, one which accounts for the average modification of the electron momentum inside the nucleus, 
and one which would be utilized for the calculation of the average focusing of the electron wave in the 
nuclear volume. 

For quasielastic (e, e') scattering a comparison of EMA calculations with numerical results from 
the 'exact' DWBA calculation |29l I3U| seems to indicate a failure of the EMA Also the improved 
approximation including a first order correction from eq. ((31 is of limited validity. However, we find 
that the EMA with an effective potential V ~ (0.75...0.8)Vo for heavier nuclei represents a viable 
method for the analysis of Coulomb distortion effects in inclusive quasielastic electron scattering, if 
the initial and final energy of the electrons and the momentum transfer are sufficiently large. 
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Figure 1: Amplitude (^7°^) 1 / 2 
of an electron wave incident on a 208 Pb nucleus with an energy of 100 MeV. The focusing varies 
strongly inside the nucleus, which has a radius of the order of 7 fm. 



2 Exact numerical calculations 

2.1 Properties of exact wave functions 

The nuclei of 40 Ca and 208 Pb were chosen for our calculations as typical examples for medium and 
highly charged nuclei. The Dirac equation was solved by using a partial wave expansion, which is 
discussed in detail in the appendix. The radial integration of the radial functions was performed by 
the method presented in [32]. However, we did not neglect the electron mass in our calculations, 
although mass effects are quite small in our case. 

We present results for an incoming electron with spin in direction of the electron momentum 
scattered off the fixed electrostatic potential of the nucleus. Considering different spin or final state 
waves with incoming spherical wave would lead basically to the same conclusions. 

The charge distribution of the 208 Pb nucleus was modeled by a Woods-Saxon distribution 

with r]/ 2 = 6.6 fm and diffusivity a = 0.545 fm, compatible with an rms charge radius of 5.5 fm and a 
central Coulomb potential of Vq = —25.7 MeV, whereas for the 40 Ca nucleus a three-parameter Fermi 
form was used 

l+uj(r/r 1/2 ) 2 

p(r) = po—, (8) 

with ri/2 = 3.766 fm, diffusivity a = 0.586 fm and oj = —0.161, compatible with an rms charge radius 
of 3.48 fm and a central Coulomb potential of V = -10.4 MeV [331133] . 

Fig. ^ shows the focusing (^7°'0) 1 / 2 of an electron wave incident on a 208 Pb nucleus with an 
energy of 100 MeV. The wave is normalized such that the density Tp^/°i/j approaches the value 1 in the 
asymptotic region. The focusing is smaller in the upstream side of the nucleus and grows larger in 
the downstream side. At the same time, there is a strong decrease of the charge density in transverse 
direction to the electron momentum. 
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Figure 2: Focusing factor of the electron wave incident on a Pb nucleus as a function of the asymp- 
totic electron momentum. 

As a first step, we checked the focusing factor in the center (r = 0) of the 208 Pb nucleus. The 
exact central focusing factor and the value typically used in the EMA match extremely well already 
at relatively low energies above 40 MeV, as shown in Fig. |2J For an electron energy of 100 MeV 
the EMA focusing factor is given by (100 + 25.7)/100 = 1.257, as shown in Fig. ©, and the exact 
value deviates less than half a percent from this approximate result. This positive result turned out to 
be generally valid for 'well-behaved' types of potentials like, e.g., Gaussian potentials with depth and 
spatial extension comparable to depth and extension of the nuclear potential. But a closer look at the 
electron-wave amplitude reveals that the amplitude varies strongly inside the nuclear volume and the 
average focusing deviates from the value calculated from the electrostatic potential Vq in the center of 
the nucleus. 

Fig. El shows the decay of the wave amplitude on an axis perpendicular to the electron momentum 
which goes through the nuclear center for an electron incident on 208 Pb with an energy of 400 MeV. 
The central focusing factor ~ 425.7/400 ~ 1.0642 decreases to 1.0187 at the outer edge of the nucleus 
at a transverse distance of 8 fm to the center. 

Plotting the focusing factor along the axis through the nuclear center parallel to the electron 
momentum shows a strong increase on the downstream side of the nucleus (Fig. 

One is therefore naturally lead to the idea to calculate an averaged focusing factor / defined by 



where p(r) is the nuclear matter density distribution which can be well approximated by the charge 
density profile of the nucleus for sufficiently large mass numbers A > 20. For a typical electron energy 
of 400 MeV, one obtains / = 1.050 for 208 Pb, corresponding to an effective potential value of —20.07 
MeV, in contrast to the often used potential depth Vq = —25.7 MeV. The same calculation for 40 Ca 
leads to an effective potential of —7.76 MeV, compared to Vq = —10.4 MeV. 
Defining an effective potential value V by 



P 



J d 3 rtp T (r)'j 'ijj T (f)p(r) 



(9) 



J d 3 rp(r) 



V = 



J ^rjj} T (r)j tp T (r)p(r)V(r) 



(10) 
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Figure 3: Focusing factor of the electron wave incident on 208 Pb with electron energy 400 MeV, on a 
straight line through the nuclear center transverse to the electron momentum. 
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Figure 4: Focusing factor of the electron wave incident on 208 Pb with electron energy 400 MeV, along 
a straight line through the nuclear center parallel to the electron momentum. 
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which is a measure for the average (semiclassical) electron momentum inside the nuclear medium, 
leads to the very similar results V = 20.12 MeV for 208 Pb and V = 7.78 MeV for 40 Ca. The difference 
between the effective momenta for the ingoing and outcoming electron calculated from this effective 
potential value can also be viewed as the effective momentum transferred by the electron to the nucleon 
in a quasielastic knockout process. 

These observations are a strong argument that a modification of the EM A with an effective mo- 
mentum corresponding to an effective potential ~ (0.75... 0.8) Vo would provide reasonable results when 
experimental data are corrected due to Coulomb distortion effects. An effective potential value of 
(18.7 ± 1.5) MeV for 208 Pb was extracted by Gueye et al. [22] by comparing data from quasielas- 
tic positron scattering to quasielastic electron scattering data taken by Zghiche et al. 36 . Exact 
calculations which include the Coulomb distortions of positrons will be presented in a forthcoming 
paper. 

It is interesting to note that also for lighter nuclei like 40 Ca a similar effective potential value should 
be used like in the case of heavy nuclei. But in such cases, Coulomb distortions are usually of minor 
importance for the data analysis, and the choice of the effective potential value that is used in the 
EMA analysis plays a minor role. 

We further mention that the average potential inside a homogeneously charged sphere is given by 
4Vo/5, where Vq = —3aZ/2R is the value of the potential in the center of the sphere. For positrons, 
we found that the same effective potential (with opposite sign) can be used, since the absolute values 
of the effective potentials differ by less than 0.05 MeV for an energy range of several hundred MeV. 

The first order term in eq. (0) fails to provide a satisfactory picture of the focusing inside the 
nucleus. It does not reproduce the strong increase of the focusing on the downstream side of nucleus, 
and the decrease of the focusing in transverse direction is also not contained. On the contrary, the 
dominant imaginary term ib(k' xr) 2 causes an increase of the modulus of in transverse direction. 

The first order term g^ accounts for the deformation of the wave front near the nuclear center, but 
higher order terms are needed in order to describe correctly the amplitude of the distorted electron 
wave inside the nucleus. Calculations with a phenomenological expression for the second order term 
g( 2 ) have been presented in JH]- The dash-dotted lines in Figs. El and 0] show the focusing that would 
be obtained from the first order expression (k'/k)\l + g^'\ for a homogeneously charged sphere with a 
radius R = 7.1 fm and Z = 82, which is a good approximation for a 208 Pb nucleus. 

It is solely the air 2 -term which accounts for a decrease of the focusing in transverse direction, 
but even if one neglects the other terms in g^ which cause an increase of the focusing in transverse 
direction, the air 2 -term leads to a negligible effect compared to the actual transverse decrease of the 
Coulomb distortion. Therefore the assumption was made in |25U26| that the focusing is nearly constant 
in transverse direction; the results presented there should be corrected for the overestimated focusing. 
However, the results in [25] are in good agreement with the exact calculations presented by Kim et al. 
[H], where exact electron wave functions were used. 

A better approach than given by expansion (^Q) to take the local change in the momentum of the 
incoming particle into account is to modify the plane wave describing the initial state of the particle 
by the so-called eikonal phase Xi (r) ( see 1^3 EZ| and references therein) 

e ikif _^ e ikir+iXi(r) ^ 

where 

2 

Xi(r) = ~ J V{f+ kis)ds = - J V(x, y, z')dz' (12) 

— oo — oo 

if we set ki = k\e z . In analogy to eq. (JHJ), the z-component of the momentum then becomes in eikonal 
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approximation 

kze ikiz+i Xi = _ idze ikiz+i X i = (£i _ v)e ik * z+iXi . (13) 
The final state wave function is constructed analogously 

e ik f f-iXf(r) } (14) 

where 

oo 

Xf (r) = - [ V(f+ k f s')ds' . (15) 



o 

In order to check the quality of this approximation, we calculated the phase of the first (large) com- 
ponent ipl/ 2 from the exact electron spinor along the z-axis, and extracted the quantity xf x by setting 

e^ z+i ^=^l /2 (z)/\^ /2 (z)\. (16) 

If the eikonal approximation (fT2"|) were exact, then the derivative ^Xi X ( z ) would be equal to the 
negative value of the electrostatic potential 

^-xT(z) = -V(z). (17) 

In Fig. © a comparison of -j£Xi X ( z ) to the absolute value of the electrostatic potential of 208 Pb is 
shown. For energies above 300 MeV, the phase of the electron wave function is very well described by 
eq. (|12j) inside the nuclear volume. 

The validity of our calculations was verified by reinserting the electron wave functions into the Dirac 
equation. This way, the electrostatic potential can be reproduced and compared to the initial potential. 
Additionally, we checked current conservation, the asymptotic behavior of the wave functions far away 
from the nucleus and the behavior for the limit Z — > 0. Finally, the calculations were also performed 
for the potential of a homogeneously charged sphere. In this case, analytic expressions are available for 
the radial wave functions occurring in the partial wave expansion |38| . which were in perfect agreement 
with the results obtained via radial integration. 

2.2 EM A from DWBA 

We establish now a connection between the DWBA results and the EMA. The DWBA transition 
amplitude for high momentum transfer in inelastic electron scattering can be written for one-photon 
exchange as 

T if = d 3 r e $r N \p e (r e )pif(r N ) - j e (r e )Jif{f N )\-^ — (18) 

J L 1 i \r e — rN\ 

with p e , pif, j e , Jif being the charge and current densities of the electron and the nucleus, respectively, 
and uj is the energy loss of the electron. The double volume integral presents a clear numerical 
disadvantage of this expression. According to Knoll [SHI, one may introduce the scalar operator 

^ = ^Ef^ ± #V^ (19) 



n=0 



Q — Ul 



such that the transition amplitude can be expanded in a more convenient form (Q = q — u> 



_4vr 
- if ~ Q 2 



d s r 



p if (r)Sp e (r) - J if Sj e (r) . (20) 
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Figure 5: Comparison of (derivatives of) the phase of electron wave functions along the z— axis for 
electrons incident on a 208 Pb nucleus with different energies. 



The single integral is limited to the region of the nucleus, where the nuclear current is relevant. The 
expansion (|19|) is an asymptotic one, which means that there is an optimum number (depending on q) 
of terms that give the best approximation to the exact value. Considering terms up to second order in 
the derivatives only one obtains 



i + 



2iq V + A 4(q V) s 



Q 2 



(Q 



2\2 



e- tqr Pe 



r) + current terms 



(21) 



Note that the e~ iqi '-term in front of the electron charge or current density cancels the spatial oscillatory 
behavior of the densities in the asymptotic region far from the nucleus, such that the gradient operator 
V in eq. ()21j) probes only the distortion of the electron current due to the nuclear Coulomb potential. 
Therefore, eq. (|21|) provides a satisfactory approximation for many interesting cases when Q 2 ^> qAk 
is fulfilled. E.g., for a typical initial electron energy ej = 400 MeV used in inclusive quasielastic electron 
scattering experiments and energy transfer uj = £j — £/ = 100 MeV, the contribution to the cross section 
due to the (q V) 2 -term in eq. (|21(1 is only of the order of 2% for momentum transfer q > 350 MeV 

0113]. 

Expansion (|21j) allows for a clear comparison of the plane wave Born approximation, the EMA and 
the DWBA. The effect of the electrostatic nuclear field on the matrix element Tjj is obviously given 
by a modification of the free electron charge and current densities p e , j e via the replacements 



[p e (r)Je(r)] -> [p'e(r)J' e (r)} 



1 + 



2iq V + A 4(<f V) 2 



Q 2 



(Q 



2\2 



^ r [ Pe (r)J e (r)]. 



(22) 



Having exact wave functions at hand, it is possible to perform a numerical analysis of this modification. 
For this purpose, it is advantageous to introduce the effective four-momentum transfer squared given 



by Q' 2 = {k'i — k'j) 2 — uj 2 . Then a short calculation shows that {ei f S> 



m) 



Q 2 



K i K f 
kikf 



(23) 
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i.e., when cross sections are calculated using the EMA, the enhanced photon propagator appearing in 
the matrix element Tjj exactly cancels the focusing effect of the initial and final state wave function. A 
numerical analysis shows that this is indeed also true for the DWBA case to a high degree of accuracy. 
E.g., a calculation for relatively low electron energies 6j = 300 MeV, tt = 200 MeV and scattering 
angle $ e = 60° of 

- 2 = J d 3 r\p' e {r)\p{r) - 2 [ d 3 r\p^ ee {r)\p{r) 

1 Jd 3 rp(r) ' /free Jd 3 rp(r) 1 ' 

in analogy to eq. © leads for 208 Pb to the numerical result f' 2 /f^ ree = 0.985, i.e. the wave function 
focusing effect of ~ 17% in the electron charge density is in fact overcompensated slightly by ~ 1.5% 
compared to the free plane wave charge density due to the enhanced momentum transfer. The situation 
is similar for a backscattering angle $ e = 143°, which has also been used in the experiment |36j . 
Numerical calculations reveal that the same observation applies to the components of the current 

— 2 — * 

density. Defining focusing factors f' k for each current component j e = {j x ,jy,jz) again leads to the 
result that the effect of the S'-operator is just a cancellation of the average focusing effect in the initial 
and final state wave functions. 

Finally, we investigated the modification of the momentum transfer inside the nucleus. For this 
purpose, we defined the local momentum transfer q(r) according to 

Pe(r) = \p' e m m? (25) 

and calculated the quantities 

ql _ I d?r\q{r)\p{r) ^ = J d 3 r\q(f)\ 2 p{r) 
1 / d 3 rp(f) ' 2 J d?rp(f) 

The results for q[ and q' 2 are again compatible with an average potential of 20 MeV. Replacing p' e by 
p e in eq. (|25|) has nearly no influence on the result for the corresponding effective momenta q\ and q2, 
and again the situation is completely analogous for the current density. 



2.3 Coulomb corrections for a harmonic oscillator single particle shell model 

We point out that the discussion presented above is a most general one as long as detailed properties 
of the nuclear current are ignored. The hope is of course that the presence of, e.g., 208 nucleons 
in a 208 Pb nucleus leads to a smooth nucleon momentum distribution, such that the features of the 
individual nucleon wave functions are averaged out to some degree in the momentum regime relevant 
for our discussion. Calculations in the framework of the eikonal approximation as presented in |25j . 
but with the electron amplitudes adapted to the amplitudes from exact calculations, also agree well 
with the EMA, a fact which supports the assumption that a semiclassical description of Coulomb 
corrections can be found at high momentum transfer. Still, one cannot exclude from the first that the 
semiclassical picture of the transition amplitude which has been presented in the previous discussion 
is flawed by phase effects which enter the cross section, leading to unexpected deviations from the 
idealized EMA-behavior. 

We therefore present results for a Pb model nucleus consisting of 82 scalar protons with harmonic 
nucleon wave functions for some typical kinematical situations, which indeed show that the modified 
EMA, using an average potential value, reproduces the Coulomb effects with high accuracy for (e, e') 
scattering, if the momentum transfer and the energy of the outgoing electron are high enough. The 
smaller contribution of the neutrons was neglected in the calculations, since we will focus on general 
considerations in the following. A detailed study using relativistic nucleon wave functions including 
spin in conjunction with more realistic current models will be presented in detail in a forthcoming 
paper. 
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At the lowest energies that have been used in the positron experiment by Gueye et al. [55] . i.e., 
for an initial electron energy e~ = 224 MeV and initial positron energy ef = 262 MeV, the situation 
is a bit more involved, and we found a relatively large disagreement between the exact inclusive cross 
sections and the values calculated from the EMA with an effective potential value of 19 MeV, as will 
be demonstrated below. This observation is not astonishing, since the typical energy of the outgoing 
electrons is only of the order of 100 MeV in this case (see also Fig. 5). 

The radial wave function R n \ of a particle in the harmonic oscillator shell model potential V(r) = 
—vq + ^Mu> 2 r 2 , where M is the proton mass, is given by (n = 1, 2... and I = 0, 1...) 

^ ,(r, ^M£) iL+i SK 44 ' * = -k (27) 

such that R n i fulfills the normalization condition / °° R^ir) r 2 dr = 1 and the integral (nl\r 2 \nl) has 
the value (N = 2n + I - 2) 

(r 2 )= R 2 nl (r)rUr = x 2 (N+-), (28) 
JO 2 



with normalization constants 



2n+i+i 

Nnl = 0F(n-l)!(2n + 2/-l)!f (29) 



The Legendre polynomials are given by 



ii 



L l n-\(z) - x; ^ k *j { 2 i+2k + i)\i z * (30) 

and the energies are En = —vq + huj Q {N + 3/2). A popular phenomenological choice for the oscillator 
strength is huj ~ 41 MeV A -1 / 3 . The (n,l) = (l,5)-shell corresponds to the l/i 11 / 2 -shell in the spin- 
orbit interaction model and contains 12 protons (but 2(2/ + 1) = 22 available states, when an artificial 
spin factor of 2 is included), whereas the the (n,l) = (1,0), (2,0), (3,0), (1,1), (1,2), (2,1), (2,2), 
(1,3), and (l,4)-shell are fully occupied. We therefore assumed for the calculation of cross sections 
that the (l,5)-shell is completely filled, i.e. spherically symmetric, and adopted a weighting factor of 
12/22 leading to an rms charge radius of the Pb model nucleus of r rmg = ^/393/82xo, as can be derived 
from eq. 1)28(1 . Accordingly, x Q = 2.515 fm was used in the calculations, such that the rms charge 
radius for the nucleus adds up to 5.505 fm. 

The proton transition current density was calculated from the free form of the Klein-Gordon ex- 
pression 



Jp 



te 



2M I 



{$}f$ n lm ~ $n*mV$}} (31) 



for initial states & n i m (r) = R n i(r)Yi m (f/r) and final plane wave states $j = e f r . The transition 
charge density was calculated from the exact current conservation relation 

iojf p = Vj p , (32) 

where u = q — ej is the energy transferred to the proton as defined before. The model used for the 
proton current is simple, but it has the advantage that it allows to perform analytical calculations for 
the electron plane wave case, which permit a verification of the accuracy of the numerical calculations. 
It also captures the most important features of the (e, e') cross sections in the vicinity of the quasielastic 
peak. Furthermore, it has to be pointed out that the Ohio group used a single particle shell model, 
where the nucleon wave functions were obtained by solving the Dirac equation for each shell nucleon 
with phenomenological S-V-potentials. But a comparison of measured data and calculations shows 
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large discrepancies especially at higher energy transfer to > 150 MeV, where correlation effects and 
pion production become increasingly important. Therefore also a single particle shell model with 'exact' 
nucleon wave functions cannot be considered as an 'exact' model for inclusive quasielastic scattering. 
Hence, the optimal strategy is to find a general, model independent method which makes it possible to 
include the Coulomb distortion effect in the analysis of experimental data. 

For lower electron energies, where the EMA cross sections (cr EMA ) start to deviate significantly 
from the exact ones (a cc ), it is still possible to match the EMA and exact cross sections by using 
an effective potential value Vfn that differs from the commonly used 18 — 20 MeV for Pb. E.g., for 

= 224 MeV, o cc and o EMA are very close for Vfu ~ 27 MeV in our model for an energy transfer 
larger than 100 MeV. This effective value is very stable under distortions of the scalar proton current 
model; using a 'wrong' value for co in eq. (j32j) which is 20 MeV larger or smaller than — e/ or 
including a potential term of similar order in eq. (|31j) changes the cross sections, but the effective 
potential value does not change significantly. This is a positive result, since it shows that our general 
considerations are not too strongly model dependent. Binding energies, nuclear potentials and the 
corresponding exact nucleon wave functions are relevant for an accurate modelling of amplitude, width 
and position of the quasielastic peak, but the impact of the Coulomb distortion of electrons on the 
quasielastic cross section expressed by a fitted effective potential Vfu can also be studied to some level 
using a simplified model with a quasielastic peak that incorporates approximately the properties of the 
true quasielastic cross section. Since binding energies are indeed neglected in our theoretical model, a 
shift of the quasielastic peak by about 25 MeV to lower energy transfer is observed. 

It is very instructive to investigate the total response function S tot , which is defined in plane wave 
Born approximation by 

<h(7 ^"" vmouxS^AWI^Q), (33) 



dn f de f PWBA 
via the well-known Mott cross section 

OMott = 4a 2 cos 2 (G/2)ef /Q 4 (34) 

not only for the total cross section, but also for individual filled shells. According to eq. (|23|) . the 
Mott cross section remains unchanged when it gets multiplied by the EMA focusing factors and the 
momentum transfer Q 4 is replaced by its corresponding effective value. Therefore, if the EMA is a 
good approximation, the cross section can be written (see also |35| 1 

cPa^r, d 2 a 



CC EMA _ atot 



dVLfdef d^lfdef 



c-Mott*Sf WBA {\q efS \ } u:,Q). (35) 



Table 1 shows ratios of total responses, which have been obtained by dividing the same Mott cross 
section out of the inclusive Coulomb corrected, EMA-, and PWBA-cross sections restricted to single 
closed shells. The first two columns show a case where both the initial and final electron energy is 
larger than 300 MeV, the other four columns show two kinematical settings where both the initial and 
final electron (positron) energy is smaller than 300 MeV. In all three cases, an effective potential value 
of 19 MeV was used. 

It is an interesting point is that even closed shells show an EMA-like behavior at higher electron 
energy. This is indeed not the case for single states, which are not spherically symmetric. Summing 
over all shells, one obtains a ratio S t °^/S t °^ I . which is practically one for an effective potential value 
V = 19 MeV in the case where = 485 MeV. Note that the different shells contribute differently to 
the total cross section, according to the number of protons contained in them and their momentum 
distribution. E.g., the (l,0)-shell, which show the most irregular behavior for the kinematics shown 
in Table 1 due to its narrow momentum distribution, contributes only marginally to the total cross 
section because the shell contains only 2 protons. 
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(n, I) 


qtot,nl 1 qtot,nl 

^CC ' EMA 


qtot,nl I qtot,nl 

CC ' PWBA 


qtot,nl I qtot,nl 
^CC ' EMA 


qtot,nl j qtot,nl 

CC ' PWBA 


qtot,nl j qtot,nl 

CC ' EMA 


Qtot,nl J qtot,nl 
CC 1 PWBA 


fi n\ 
(1>0J 


1 1 m 
1.1U1 


1 Oil 


1.2/5 


3.5/4 


0. / /6 


U.425 


(2)°) 


1.UU4 


U.9 1 1 


U.99o 


0.9/2 


1.056 


U.D /9 


(3,0) 


l.Ulo 


1.U40 


1 1 1^ 
1.110 


1.2 / 2 


1.033 


U.8 /U 


(1,1) 


1.024 


1.332 


1.140 


2.041 


0.896 


0.758 


\ L > Z ) 






l.UDo 






u.yoz 


(2,1) 


1.013 


1.268 


1.050 


1.377 


1.025 


0.718 


(2,2) 


0.995 


1.271 


1.031 


1.526 


1.028 


0.939 


(1,3) 


0.982 


0.989 


1.054 


1.132 


1.003 


1.087 


(1,4) 


0.990 


0.956 


1.072 


1.043 


1.011 


1.115 


(1,5) 


1.005 


0.955 


1.085 


1.011 


1.006 


1.117 


£ 

(n,0' 


0.999 


1.072 


1.073 


1.244 


0.989 


0.956 



Table 1: Ratio of the total response for different filled shells, calculated in plane wave approximation 
(S t °J VBA ) and EMA with an effective potential 19 MeV (S^ f J, with the exact response S^. First 
two columns: Initial electron energy = 485 MeV, uj = 140 MeV, ?? e = 60°. Columns 3 & 4: Initial 
electron energy = 224 MeV, uj = 100 MeV, $ e = 143°. Columns 5 & 6: Initial positron energy 
€i = 262 MeV, uj = 100 MeV, ~Q P = 143°. The bottom line shows the ratios of the total response 
functions for the full cross sections, i.e. after summing over all shells with a weighting factor 6/11 for 
the (l,5)-shell. Note that these values are not averages of the data in the column above, since the 
contribution of each shell to the total cross section is different. 

In the case given in Table 1 with an initial electron energy e~ = 224 MeV, the exact cross section 
is 7.3% larger than the EMA cross section and 24.4% larger than the plane wave cross section. By 
naive interpolation, one can infer that the exact and the EMA cross section should be identical for an 
effective potential value which is larger than 19 MeV. This indeed the case for an approximate value 
Vfu — 27 MeV, and similarly in the positron case for Vfu — 25 MeV. The maximum values of the 
total response for the e~ = 224 MeV and ef = 262 MeV case differ by about 10%, i.e., one observes a 
relatively large deviation from the EMA prediction that the two responses agree. However, for higher 
energies used in the positron experiment of Gueye et al., an effective potential value of 18.9 ± 1.5 MeV 
is compatible with our model calculations. 

The accuracy concerning the calculation of cross sections and accordingly the values given in Table 
1 is limited to about 1% due to the truncation of the Knoll expansion eq. ()21|) and the finite resolution 
of the grid that has been used for the modelling of the nucleus. The numerical evaluation of transition 
amplitudes was performed by putting the nucleus on a three dimensional cubic grid with a side length 
of 30 fm and a grid spacing of (30/140) fm, and convergence was checked by using different side lengths 
and grid resolutions. The accuracy of the solid angle integration of the (e, e'p) cross section was better 
than 0.05%. 
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3 Conclusions 



The complex behavior of the Coulomb distortion of electron waves at relatively high energies rel- 
evant for quasielastic electron-nucleus scattering experiments was studied using accurate numerical 
calculations. Naive lowest order approximations in aZ are not suitable for the analysis of Coulomb 
corrections in scattering experiments, unless they are modified in a well-controlled manner based on 
exact calculations. A Fortran 90 program is available now which can be used for accurate calculations 
of continuum electron wave functions in a central electrostatic field. The exact wave functions were 
used for a numerical study of Coulomb distortions in inclusive quasielastic electron scattering. It turns 
out that the effective momentum approximation is not reliable when the central potential value Vq of 
the electrostatic field of the nucleus is taken as basis for the EM A calculations, but a smaller average 
value V ~ (0.75...0.8)Vo leads to very good results, if the momentum transfer and the energy of the 
scattered electron large enough. An effective potential value of 19 MeV is a very good choice for 
ef > 300 MeV and Q 2 > (400 MeV) 2 , and we conjecture that at very high energies, a limiting effective 
potential value close to 20 MeV is reached. If the energy of the scattered electron becomes smaller 
than 300 MeV, the semiclassical description of the final state wave function becomes obsolete, but it is 
still possible to use an EMA-like approach for the description of the inclusive cross section by using a 
modified fitted potential value, given the condition that the initial and final energy of the electron and 
the momentum transfer are not too small. In this region, detailed calculations become necessary with 
more refined nuclear models that the one used in this work. However, for the kinematical settings that 
will be used in the future experiments at the TJNAF, our analysis shows that the EMA will provide 
a valuable strategy for the correction of data. 

It highly improbable that the match of our exact cross sections with the EMA cross sections (i.e., 
with V = 20 MeV for 208 Pb), which is better than 2% for the kinematical region e f > 300 MeV 
and Q 2 > (400 MeV) 2 , is just a pure coincidence, since both approaches are unrelated and based on 
different calculational strategies. We must therefore conclude that our findings are not compatible 
with the conclusions drawn in jSj. 
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Appendix: Calculation of continuum states of an electron in a central 
electrostatic field 

Some details concerning the description of electron continuum states in a central electrostatic field 
are given here for the reader's convenience and in order to give a fully consistent description of the 
problem. 

We first consider the solutions of the stationary Dirac equation 

[-iaV + mj3 + V{r)]ip(r) = E^(r) (36) 

for an electron with total energy E subject to the central electrostatic potential generated by a spheri- 
cally symmetric nucleus with charge number Z. We use standard Dirac and Pauli matrices [301 • Dirac 
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spinors describing states with definite angular momentum and parity can be decomposed into a radial 
and an angular part 

~ { iUr) X t K (r) ) ^ ] 
with two-component spinors Xk which are eigenstates of the spin-orbit operator 

(at + 1) X % = -«x£ (38) 

and the angular momentum operators 

P x Z = (L + s) 2 xZ = j(j + l)xt, 

L 2 X £ = l(l + l)xZ, s 2 x» = \x^ 

■M = MX£, (39) 
where k = ±1, ±2, ... is related to j and I by 

K = 1(1 + 1)- (j + ^) 2 , J = \k\-^, l = j + ^sgn( K ), (40) 

and the operators L and s are given by L = —ifx V and s = \o. The spinors can be expressed using 
Clebsch-Gordan coefficients as 

x»= E a(M-C)kb»^r c ( f )xc, 3 = V¥+i (4i) 

<=±l/2 

where XC are standard Pauli spinors, or more explicitly 



X£ = I "Sv ^ j for k < 0, = ( I for « > 0. (42) 

The radial functions / re and g K fulfill the coupled differential equations 

(43) 



±9k) = E + m-V) g K 

dr[ f K { -(E-m-V) ^ [ f K 



For an electron with energy E > in the Coulomb field of a point like charge eZ the potential V is 

V(r) = -e/r, £ = "Z, (44) 
and the continuum solutions of the Dirac equation are given by (see 0^ and references therein) 



/7k 



Re 
Im 



[(7k + i7?)e^e- ifcr F( 7K + 1 + ir/, 2 7k + 1; 2*fcr)], (45) 



where 
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e2iip= -K + ir,m/E k = ^T^. (46) 

Ik + iV 

ip is positive for aige 2lip G [0, vr] and negative for arge 2 ^ G (— vr,0]. These Coulomb wave functions 
have the asymptotic forms (r — > oo) 

1 7T 

0«(r) ~ T-coB(fcr + 77log2*T-(i + l)- + (5«) J (47) 



/ K (r) ~ -JE-JH S in(/cr + r ? log2A:r - (/ + 1)£ + £„) , (48) 
y hj -\- ui Kt z 

where 

1 —K H fa- 7T 7T 

°k = ^arg — argT(7 K + ir]) - j K ~ + (I + 1)- := 

^-argr( 7 + ^) + (/ + l- 7 )|. (49) 
In the limiting case Z —> 0, we have 



5 re (r) -> -sgn(K)ji(A:r) , / re (r) ->■ -W ^ m j v {kr), (50) 

\ Hi -\- m 

where ji(kr) = Jl^t ^+1/2 (^ r ) are spherical Bessel functions of the first kind, and V = j — ^sgn(/t) = 
1(—k). Replacing 7 K by — 7 K in eq. (|4*5|) and correspondingly in eq. leads to the irregular solutions 
g l K , f l R . Their asymptotic forms are also given by eqns. ([47148)1 . if one replaces the expression for the 
phase shift <5 k (7k) of the regular solutions by 5' K {^ K ) = 5 K (—j K ). 

The calculation of the wave functions g K and f K is most simply performed by using the real series 
expansion 

( 9 f)= (M*" 1 £ ( 12 ) (*r)» . (51) 

Inserting the series expansion ()51j) into eq. (fiU)) leads to the coupled recursion relations [3J| 

(n + l)(2 7re + n + l)ka K - n+ \ + £{E - m)a K - n 

-(j K + n+l- k)(E + m)b K , n = 0, (52) 



(n + 1)(2 7k + n + l)fcfe K;n+ i + f (E + m)6 K;n 

. (53) 

The series expansion ()51j) can be shown to converge for all values of r. 

We give here some details for the derivation of the recursion relations. From 

— g(r) = (E + m + £/r)/(r) - ^-g(r) (54) 
ar r 



one readily derives 



(7 - l^kry^k anikrT' 1 + (7 - l^kry 1 — + {kr^k £ na n (fcr) n ~ 1 

n=l ^ n=l 



(£ + m)(A;r) 7 - 1 ]T 6 n (A;r) n + ^fer) 7 "^ £ ^(Ar)" -1 + £(jfer 



-1&0 



n=0 n=l 
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y-l a O 

■JL n \KT) — T LjyKTj 

n=l 

where we have omitted the index k for notational convenience. Comparing the lowest order terms 



(k + l){kry~ l k anikr) 71 ' 1 - (k + l)^)^ 1 ^, (55) 

3 have omitted tl 
~ r 7_2 immediately leads to the starting relation 



6 = — "z — ao = o . (56) 

f « -7 

From 

A /(r ) = _(£ _ m + Z/ r )g{r) + — /(r) (57) 
ar r 



one obtains 

( 7 - l)(kry- l k bnikrY 1 - 1 + (7 - ^(fcr) 7 " 1 - + {kr^k ^ n^Ar)"" 1 



n=l n=l 



-(£ - m)^) 7 " 1 £ a n (kr) n - ^kr^k a^kr)' 11 ' 1 - ^Ar) 7-1 — 

n=0 n=l ' 



oo I 

+ (k - lXfcr) 7 - 1 *; ^ 6 n (/tr)"- 1 + (k - l)^)^ 1 -. (58) 

71=1 r 

Comparing again the lowest order terms ~ r 7 ~ 2 leads again to the starting relation (|56|), Taking into 
account the higher order terms ~ (A;r) 7_1 (/cr) n_1 in eqns. (|55I58|) leads to 

(7 — l)ka n + nka n — (E + m)b n _i — £,kb n + (k + l)&a n = (59) 

(7 — l)fcfc n + nkb n + (E — m)a n _i + £/ca n — (k — l)M> n = 0. (60) 
Replacing b n from eq. ((ST))) 

7 + K + m E + m 
b n = 1 a n -r—b n -i (61) 

in eq. (|H0|) gives 

n(27 + n)fco n + £(E - m)a n _i — (7 + n — k)(E + m)b n -i = 0, (62) 

which is equivalent to eq. (|52|). Recursion relation (|53f) is obtained analogously. 

An incident (outgoing) electron with asymptotic momentum k, energy E and polarization r is given 
by a linear combination of the 



W = E+Jl J2 e ±i5 H\l (m - r) i t I J /i)(!T T (k)mr-). (63) 

It is instructive to consider, e.g., the first component of the Dirac spinor ip T for an electron incident 
along the z-axis with spin in the same direction. Then, the term (y^ _T (k)* is only non-zero when 
\i = r = 1/2, and we have 



Y?($M = J^fl(co80), = o ttp = o) = J*+± (64) 
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Therefore, the first component of tp T is given by 



^ = v^^^ (2Z + 1)e "" <I<Z0 ^ l 4 )2 ^ (r)fl(co8,?) - (65) 



A straightforward calculation shows that 



11, -1,2 



(10 — \j-r = • (66) 

x 22 IJ 2 / 2k + 1 v ' 

Furthermore, the asymptotic behavior of g K is given by 

1 1 —K + ir/m/E ir 

9kv) ~ -r- cos(kr + i]log2kr — argT(7 + irj) H — arg 7— ). (67) 

kr 2 7 + ir\ 2 

We consider now the limit Z — > for this asymptotic expression. For Z — > we have also rj —* 0, i.e. 
argT(7 + ir/) — ► and 7 — > \k\. From 



2iu> -K + irjm/E -Kj + rjm/E ,r]K + r]m/E 

e * = ; = 5 5 h i 5 5 — (68) 

7 + zr? 7 Z + 7 + 



we see that the argument of ip approaches ir/2 for ij —* and k > 0, whereas for k < we have — > 0. 
This shows that the asymptotic behavior of the g R for Z — > is given by 

5k (^) ~ 1— cos(£x - — + sgn(«)-) = -sgn(/c) sin(/cr - — ), (69) 
kr 2 2 2 

in accordance with eq. (|5(J[) which states that the free- field solutions of g R are given by g K (r) = 
—sga(K)ji(kr). Therefore, the terms in the expansion (|H5j) for k > become (/ = k, 5 k — > tt) 



1 E + m 



J2H l e in ^ji(kr)Pi(cos^) (70) 



2E ^-i kr 

K>0 



and for k < (/ = — k — 1, S K — ► 0) 

/i: + ^ ^(/ + l)i'^(^)PKcos^), (71) 



2£ ^ v ' At" 



i.e. we obtain the partial wave expansion of a plane wave, and the normalization is such that the full 
free spinor for arbitrary momentum k and helicity £ is given by 



For the case of a realistic nuclear electrostatic potential, analytic expression for the radial functions 
are no longer available. Therefore, we calculated the radial wave functions by numerical integration 
according to the method described in appendix 3 of [32]. Outside the nuclear charge distribution 
(i.e. for r > 14 fm in our actual calculations), the electrostatic potential is a Coulomb potential, and 
therefore the radial functions G K obtained from the numerical integration can be written as a linear 
combination of regular and irregular solutions of the Dirac equation with a Coulomb potential 

G K = c K g K + d K gl (73) 
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Since the asymptotic behavior of the regular and irregular radial functions is given by 

1 7T 

9k (r) = — sin(£x + r]\og2kr - I - + S K ) (74) 

1 7T 

9 K ( r ) = -j^sm(kr + r]log2kr - I- + 5' K ) 

= S i n (kr + v \og2kr- + S K + (S' K - S K )), (75) 
kr 2 

and the asymptotic behavior of G K is described by the phase shift A K via 

G K (r) ~ A sin ( z + A K ) (76) 
with z = kr + r\ log 2kr — I ~ + 5 K , such that we obtain for A re 

A sin(z + A K ) = A sin(z) cos(A K ) + A sin(A K ) cos(z) 

= c K sin(z) + d K sin(z + 8' K - 8 K ) 
= c K sin(z) + d K sin(z) cos((^. - 5 K ) + d K sin(^ - 5 K ) cos(z), (77) 

the relations 

A cos A K = c K + (i K cos(<^. - 6 K ) , A sin(A K ) = d K sva{5' K - 5 K ) (78) 

and therefore 

sin(<5' — 5 K ) 

tan A K = V g " ; , 79 

c re /d re + cos(d' K - d K ) 

which fix uniquely the phase shift e l ^ K . The radial function G K (and the corresponding F K for the 
lower spinor components) obtained from the numerical integration procedure must subsequently be 
multiplied by a factor A -1 , where 

A 2 = (A cos A K ) 2 + (A sin A K ) 2 = (c 2 + d\ + 2c K d R cos A K ) (80) 

such that 

G K = \G H (r) ~ ^ sin(A;r + A K ) (81) 
A kr 

is correctly normalized. 

The expansion for an incoming wave scattering off a spherically symmetric nuclear charge distri- 
bution is finally given by 



* T = W £ (M - r) i t I j ^)yr r *(k)C(r), (82) 

where V'k is defined according to eq. (|57j) with the radial functions g K and f K replaced by G K and F K . 
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